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We derive the effective energy-momentum tensor for cos- 
mological perturbations and prove its gauge-invariance. The 
result is applied to study the influence of perturbations on the 
behaviour of the Friedmann background in inflationary Uni- 
verse scenarios. We found that the back reaction of cosmolog- 
ical perturbations on the background can become important 
already at energies below the self-reproduction scale. 
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Introduction 

It is well known that gravitational metric perturba- 
tions treated as propagating on a curved "background 
space-time" have an effect on the evolution of this "back- 
ground" . This is due to the nonlinearity of the Einstein 
equations. A convenient way to describe the back re- 
action of fluctuations on the background is to consider 
the "effective" energy-momentum tensor (EMT) for these 
metric perturbations. 

This problem has been studied by several authors in 
applications concerning gravity waves (see e.g. and 
references therein). One of the main puzzles in need to be 
solved is the problem of gauge invariance of the effective 
EMT. Namely, the effective EMT should be defined in a 
manner that the answer to the question "how important 
are perturbations for the evolution of a background?" 
does not depend on the choice of space-time coordinates 
(in other words, it should not depend on the gauge). 

The issue of gauge invariance becomes critical when we 
attempt to analyze how gravitational waves and scalar 
metric perturbations produced in the early Universe 
influence the evolution of the background Friedmann- 
Robertson- Walker (FRW) Universe. The procedure sug- 
gested by Isaacson Q defines a gauge-invariant EMT for 
small-wavelength, high-frequency perturbations, and is 
not applicable in our case for the following reason. In or- 
der to get the invariant EMT following this prescription, 
one should average terms in the Einstein equations which 
are quadratic in the perturbations over time intervals big- 
ger than the typical inverse frequency of perturbations. 
Obviously, it is assumed that the time scale character- 
izing the background is much bigger than the period of 
the perturbations. Since in the early Universe inhomo- 
geneities with scales bigger than the horizon scale are 



frozen, it means that their typical period is much bigger 
than the cosmic time scale and the procedure cannot be 
used. 

In this Letter we consider perturbations about a FRW 
manifold and show how to define a gauge invariant EMT 
for metric perturbations which involves only spatial av- 
eraging on a hypersurface of constant time. This allows 
us to formulate the problem of back reaction of perturba- 
tions on the evolution of the background FRW Universe 
in a coordinate-independent manner at every moment in 
time. 

We apply our framework to a chaotic inflationary 
model. Given the spectrum of linear cosmological pertur- 
bations generated during inflation, we evaluate their ef- 
fective EMT and find that back reaction becomes impor- 
tant already at energy scales lower than those at which 
the stochastic driving terms dominate. This may have 
important consequences for the dynamics of chaotic in- 
flationary models. 

There has been recent work on the back reaction of 
density inhomogeneities in cosmology. Futamase |^ con- 
sidered the problem of back reaction in harmonic gauge. 
Seljak and Hui reconsidered this issue using a differ- 
ent gauge but obtained differing results, thus highlight- 
ing the need for a gauge-independent analysis. A similar 
problem was also addressed by Buchert and Ehlers in the 
context of Newtonian cosmology . 

This Letter is organized as follows: In Section 2 we 
formulate some useful properties of the diffeomorphism 
transformations. The back reaction problem is set up in 
Section 3, where we show how to formulate it in terms 
of gauge-invariant quantities only. Section 4 contains an 
application of our results to study the back reaction prob- 
lem in the chaotic inflationary scenario. 

Diffeomorphism transformations 

The gauge group of General Relativity is the group 
of diffeomorphisms. To define it we consider a smooth 
vector field on the space-time manifold Ai . The set of 
parametrized integral curves of ^" are given by solutions 
of the differential equations 



(1) 



(A being an affine parameter) with initial conditions 
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X"(A = 0) = for every a;". This induces a coordi- 
nate transformation on Ai (see also M): 



x" — > = x"(A = 1) = e^'^if^a;" 



(2) 



where ^ should be considered small if we want to use a 
perturbative expansion in (^). 

Now let us take two different points P and P of the 
manifold A4 having the same coordinate values Xq in the 
two distinct coordinate frames x and i, that 
and x'p — Xq . We want to express the value of an 

arbitrary tensor field Qp at point P in the coordinate 
system x in terms of Qp and its derivatives at point P in 
the coordinate system x. The answer is well known and 
is given by the Lie derivative: 



Q(xo) = (e-^«Q)(a;o) 



This Lie operator obeys an important property, which 
we exemplify below in the case of the Eintein tensor G. 
We can express G as a function of the metric and its 
derivatives: 



G{X) EE G 



d_ 

dx 



9{x) 



Since the diffeomorphism transformation 
feet the derivatives one can write 



(4) 



does not ef- 



(e-^«G)(a;) = G 



d_ 

dx 



(5) 



Regarding G{x) as a Junctional of the metric we can ex- 
pand (|^) in terms of functional derivatives and obtain for 
example the following property of the Lie derivative: 



C,G{x) = j d^x'^^^C,g{x') 



(6) 



where 5G{x) / 5g{x') is the functional derivative of the 
Einstein tensor with respect to the metric. Formulas 
similar to (||) are true also for the EMT and in fact for 
arbitrary tensor fields which can be considered as local 
functionals of other tensor fields and their derivatives. 

Back Reaction and Gauge Invariance 

We consider a FRW Universe with small perturbations. 
This means one can find a coordinate system {t^x"^) in 
which the metric {gp,u) and matter fields denoted for 
brevity by the collective variable g° = (g^^,if), can be 
written as 



q'^it,x^) ^q^it) + dq'^it,x') 



(7) 



where qQ{t) depends only on the time variable and 
-C |(7o|. It is also assumed that the spatial average 
of (5g° over hypersurfaces t =const with respect to the 
induced "homogeneous" part of the 3-metric vanishes. 
The Einstein equations 



Guv - SttTm^ := Iluiy = 



(8) 



can be expanded in a functional power series in Sq"" about 
the background (t) if we treat G^i, and T^i/ as function- 
als of 9" , namely 
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n(g^) + naSq'' + -UabSq'^Sq" + 0{Sq, 







(9) 



(omitting tensor indices). From now on we adopt De- 
Witt's condensed notation i.e. assume that con- 
tinuous variables (i, a:*) are included with the field in- 
dices a, b so that, for instance, g° = q°'{t',x'^) and 
(3) n = 5Il/5q"'\^^ etc. In addition, the summation over 
repeated indices is understood to include integration over 
time and/or space. 

To lowest order, the background go(*) s-'^d the pertur- 
bations Sq"" satisfy, respectively, the equations 



n(gg) = and n,a^g" = . 



(10) 



However, it is clear from (|^) that to next order in 5q 
the perturbations also contribute to the evolution of the 
background homogeneous mode of the metric and mat- 
ter fields gg. To see this, we take the average of (^ over 
a t =const hypersurface, and obtain the following "cor- 
rected" equations for the evolution of the background: 



n(gS)= -^{'il,ab5q''6q^) 



(11) 



where brackets () denote spatial averaging. At first 
glance, it seems natural to identify the quantity on the 
right hand side of Eq. (|ll|) with the effective EMT of 
perturbations which describes the back reaction of per- 
turbations on the homogeneous background. However, 
this expression is not invariant with respect to diffeomor- 
phism transformations and, for instance, does not vanish 
for "metric perturbations" induced in Minkowski space- 
time by a coordinate transformation. 

Thus it is clear that if we want to clarify how important 
physical perturbations are for the background evolution 
we need a diffeomorphism independent (gauge invariant) 
measure characterizing the strength of perturbations. 

The coordinate transformations (H) induce diffeomor- 
phism transformations (^) on 5q which, in linear order, 
take the form 



Sq- ^ = Sq"" - C^q^ 



(12) 



where (^) — 0. To second order, the background variables 
are not gauge invariant either but change as 
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1o 



(13) 



(21) 



where 



Let us write the metric for a perturbed flat FRW Uni- 
verse 



N^{t){l + 2(f))dt^ ~ 2a2(t)(B , - S^)dx'dt (14) 



a2(t)[(l-2^)5. 



2E, 



F,_ 



hij]dx^dx-' 



where the 3-scalars 0, B, ■(/;, E characterize scalar pertur- 
bations, Si and Fi are transverse 3- vectors and hij (grav- 
ity waves) is a traceless transverse 3-tensor . 

Under a gauge transformation (|l^), the quantity = 

^^{B - E),-E^, - Fj], with a "dot" denoting time 



derivative, changes as 



(15) 



This quantity will be treated formally as a 4-vector in 
Lie derivatives below. Using one can form gauge 
invariant variables characterizing both background and 
linear perturbations: Q = e'~^ q, that is 



<5g- = <5q- + /:^g- 



(16) 



and 



Idtt 



(22) 



can be interpreted as the gauge invariant effective EMT 
for perturbations. Therefore if we want to find out if the 
back reaction of perturbations is important we should 
compare t^^{5Q) of perturbations with T^;y((5o)- Note 
that none of the terms in Equation (|2^) depends on the 
specific coordinate system used to evaluate them. 

To conclude this section, we will derive the effective 
EMT for scalar cosmological perturbations about a spa- 
tially flat FRW Universe. Since the results do not depend 
on the gauge, we can calculate the EMT using longitudi- 
nal gauge [p|, in which 



(1 + 2cj))df - a^{t){l - 2il:)5^jdx'dx^ 



(23) 



and the matter perturbation (taking matter to be a scalar 
field) is 5Lp. For many types of matter (scalar fields in- 
cluded) Tij is diagonal in linear order in 5q, which implies 
that = By evaluating the functional derivatives 

in (prj) (see also [|ll|) one can derive the following expres- 



sion for T, 



Qg - go" + {Cx5q^) + - {C\q^^) 



(17) 



It is easy to verify that the SQ'^ correspond to the set of 
Bardeen's gauge invariant variables ||]. The actually 
change under diffcomorphism transformations as 



(18) 



where X] is the commutator of the vectors ^ and X. 
For uncorrelated ^ and X we have ([^, X]) =0, and there- 
fore the last term in ( |l^ ) vanishes (see Ref. |l^ for a 
detailed discussion of this term). 

Our goal is to rewrite Equation (|l^) in terms of quan- 
tities which are gauge invariant up to second order in 
perturbations. It is easy to see from Identity (^) that if 
Einstein's equations are valid for the set of variables g, 
then 



e^-^n((7) = n(e^^(7) = n(Q) = 0. 



(19) 



Expanding (19) to second order in 5Q and taking the 
spatial average of the result yields 



n(go) = -i(n,afc<5g",5g^). 



(20) 



which is the desired gauge invariant form of the back- 
reaction equation. Note that in deriving (j2^) we made 
use of the equations of motion for q. Finally, Equation 
(pO|) can be written as 



Too = -J- 



-12H 



-3((0)2)+9a-2((V0)2) 



+ \Vvvivo){Sv^) + 2F^((^o)(0M 



(24) 



Tij — Qj Oij 



1 r 
8^ 



{24:H^ + 16H){(f>^) + 24:H 
4 



+ ((^)^)+4(00)--a-2((V0)2) 



- il^v^((^o)(V) +2F^(<^o)('/>'5¥') 



(25) 



where H — a /a is the Hubble parameter and tq^ = = 
(i^j). 

Back Reaction in Stochastic Inflation 

As an application of the formalism developed in the 
previous sections, we will evaluate the order of magnitude 
of back reaction effects in the chaotic inflationary sce- 
nario [ p^p^ , for simplicity taking a massive scalar field 
as the inflaton. In this model, quantum fluctuations of 
the scalar field certainly dominate the dynamics of the 
background when the field is above the self-reproduction 
scale ipsf ~ m^^/^ (in Planck units), and space on scales 



3 



of the particle horizon is completely inhomogeneous, con- 
sisting of many bubble Universes. It is usually supposed 
that in spatial regions where the scalar field at some point 
drops below ^sf, the evolution proceeds classically and 
the metric fluctuations generated are not very important 
for the evolution of the homogeneous background. We 
will show below that this is not really the case. 

In a chaotic inflationary universe scenario, linear per- 
turbations on a flxed comoving scale k are completely 
specified by the function (for a review, see Q). This is 
due to the fact that ijj = (p and that the metric and mat- 
ter perturbation variables (p and 5(p are anti-correlated 
for ka <C H, i.e. Sipj. ~ — (^o^fc- Hence, all terms in the 
effective energy-momentum tensor r^j, can be expressed 
through the various correlators of (pk ■ The amplitudes of 
(j)k are known from the theory of linear cosmological per- 
turbations. Using the results for (j)k valid during inflation 
M we obtain for instance the regularized correlator 



dk 
1 

to2 dk 
3271 V^(i) Jk, ~k 



(26) 



In 



H{t)a{t) 



2'^o(ii 



where t denotes physical time, ti is the time when infla- 
tion started and the inflaton potential is V = l/2m?tf'^. 
The IR and UV physical cut-offs ki and kt are given, 
respectively, by the scale of the largest wavelength per- 
turbation (created when inflation started at time ti), i.e. 
ki — H{ti)a{ti), and by the scale kt = H{t)a{t) of the 
shortest classical perturbation, which is just the scale of 
the Hubble distance. 

It can be checked that the main contribution to the 
EMT of cosmological perturbations comes from terms 
proportional to the above correlator. Therefore one finds 
that at the end of inflation (when tpQ ^ 1) the energy 
density of perturbations is about 



(27) 



Comparing the above result ( ^ ) with the background 
energy density at the same moment of time, we conclude 
that if at the beginning of infiation 



ipoiU) > (fibr ^ m 



-1/3 



(28) 



then back reaction becomes important before the end of 
infiation {ipo ~ 1). 

It is important to note that ipbr is smaller than the 
value tpsf ~ mT^/"^ when stochastic source terms from 
quantum fluctuations start to dominate. A more detailed 
discussion of back reaction will be the subject of a forth- 
coming publication 
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